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Abstract 

The new formalism for description of (un)polarized semi-inclusive DIS at mod- 
erate energies is developed. Hadron production is modeled as a product of distribu- 
tion functions and hadronization functions (HFs) weighted by the hard scattering 
cross sections, as it is done in the LEPTO event generator. The correct treatment 
of polarization effects shows that description of semi-inclusive DIS within this for- 
malism includes a new nonperturbative input - 'polarized HFs. It is shown that 
this approach does not correspond to the commonly adopted one with the indepen- 
dent quark fragmentation. The purity method used by the HERMES collaboration 
mixes up the two approaches and ignores the contributions from polarized HFs. 
This method cannot be considered a precise tool for the extraction of polarized 
quark distributions from measured SIDIS asymmetries. 

1 Introduction 

Deep inelastic scattering (DIS) is one of the main sources of our knowledge of the nucleon 
structure. More information about both the nucleon structure and a hadron production 
mechanism one can obtain by studying the semi-inclusive DIS (SIDIS). 

It is evident that the theoretical description of SIDIS is much more complicated than 
that of DIS owing to our poor knowledge of the nonperturbative hadronization mech- 
anism. Traditionally, one distinguishes two regions for hadron production: the current 
fragmentation region, xp > and the target fragmentation region, xp < 1 . The common 
assumption is that hadrons in the current fragmentation region with z > 0.2 are produced 

1 We use the standard SIDIS notations and variables like in pQ. 
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in the independent quark fragmentation. Then, in the LO approximation of perturbative 
QCD the SIDIS cross section for unpolarized target is given as 



a h (x, z, Q 2 ) cc (1 + (1 - y f) £ e\ q(x, Q 2 ) D h q (z, Q 2 ) (1) 

i 

and for polarized beam and target 

Aa h (x, z, Q 2 ) oc (1 - (1 - y) 2 ) £ e\ Aq(x, Q 2 ) D h q (z, Q 2 ) . (2) 

g 

The virtual photon asymmetry for hadron h production can be expressed as 

This equation can be rewritten as follows: 

Ai(x, z, Q 2 ) = Vi(x, z, Q 2 ) , (4) 

where the quark polarizations (Aq/q) are factored out and the purities, V q , are defined 

as 

e 2 a q(x, Q 2 ) DNz, Q 2 ) 

V h (x z O 2 ) = q 9 (5) 

M ' M ] E q/ e 2 q ,qi(x,Qi)D^z,Q*)- {b) 

Recently, the important issue of the extraction of polarized quark distributions was 
again addressed by the HERMES collaboration pQ . They have used the above LO descrip- 
tion of SIDIS and calculated purities using the Monte Carlo unpolarized event generator 
LEPTO Then, using the asymmetries, measured for different hadrons the helicity dis- 
tributions Aq(x) were extracted by solving Eq. (jlj). 

The main assumption of this method is that the hadronization mechanism in LEPTO 
is the same as in naive picture of SIDIS where all hadrons in the current fragmentation 
region with z > 0.2 are produced in independent quark fragmentation and there are no 
additional terms in both the numerator and the denominator of Eq. (jSJ). This assumption 
is based on factorization theorem of QCD which is proven as an asymptotic statement for 
very high momentum transfers and hence very high energies. 

An alternative approach is adopted in LEPTO j2j event generator. The hadronization 
mechanism of this generator is based on LUND string fragmentation model implemented 
in the JETSET program 0. In this model the QCD confinement, the quantum numbers and 
energy-momentum conservation are taken into account. As a consequence at moderate 
beam energies when the final hadronic system has a limited invariant mass, ~ 3 — 5 GeV, 
one cannot neglect the influence of the target remnant state on distributions of hadrons 
produced in the current fragmentation region. 

In Sec. El of this paper it will be demonstrated that the properties of the " quark 
fragmentation functions" extracted from generated LEPTO samples are in contradiction 
with generally accepted properties of independent quark fragmentation. The reason for 
this discrepancy is indicated. The generalization of the parton model expression for 
polarized SIDIS is given in Sec. El Finally, in Sec. 0] some discussion and conclusions are 
presented. 
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2 LEPTO and Fragmentation Functions 



In the standard picture of SIDIS (see Eq. the quark fragmentation functions by 
definition depend on the type of hadron, quark flavor and fraction of quark energy carried 
by hadron, z, (there is also a weak dependence on Q 2 due to perturbative QCD effects), 
and are independent of a) the Bjorken variable x and b) the target type. These properties 
related to universality of fragmentation functions are essential - they indicate that one is 
dealing with independent quark fragmentation and that there is no influence of the target 
remnant on hadron production in the current fragmentation region. 

These fragmentation functions are not well known for different hadron and quark types 
and the LEPTO event generator is used by HERMES collaboration P to calculate purities. 
However, hadronization in this generator is based on the LUND string fragmentation 
model and one has first to check if the issues a) and b) are satisfied in this approach. To 
this end samples of SIDIS events were generated for HERMES experimental conditions 
using the settings of LEPTO as in [T] (see also reference [64] of [I]). The option LST(8)=0 
of LEPTO was used since it corresponds to the LO approximation of SIDIS 2 . The following 
cuts are used: Q 2 > IGeV 2 , W 2 > lOGeV 2 , y < 0.85, 0.023 < x < 0.6, E' > 3.5 GeV, 
z > 0.2 and x F > 0.1. 

In Fig. [I]t he quark "fragmentation functions" to tt + , extracted from the sample, gener- 
ated for the HERMES kinematics on a proton target, are presented as a function of z. The 
available range of the Bjorken variable, x, was divided into two equally populated inter- 
vals: 1) (x < 0.094) and 2) (x > 0.094). As one can see from Fig. [H the extracted "quark 
fragmentation functions" happen to be strongly dependent on the Bjorken x variable in 
striking contradiction with the property a) mentioned above. Note, that this cannot be 
attributed to the (weak) Q 2 -dependence of fragmentation functions. To demonstrate this 
the LO fragmentation functions from ^Ej, which includes the QCD evolution, are also 
presented in the same figure for mean values of Q 2 corresponding to Bjorken variable 
intervals 1) (Q 2 = 1.5(GeV/c) 2 ) and 2) (Q 2 = 3A(GeV/c) 2 ). 

The " quark fragmentation functions" obtained from the samples generated for proton 
and neutron targets with the cut x > 0.1 are presented in Fig. |21 We see a dependence 
on the target type which contradicts property b) of fragmentation functions. 

Such behavior of "fragmentation functions" extracted from the generated samples is 
also observed for the production of other types of light meson like tt~ , K + , K~ etc. 

At this point one can conclude that the "quark fragmentation functions" extracted 
from the samples generated for HERMES kinematical conditions do not correspond to 
the commonly used notion of fragmentation function. 

Let us recall that hadronization in the LEPTO event generator is based on string frag- 
mentation and as it is stressed in j^j: "the primary hadrons produced in string fragmen- 
tation come from the string as a whole, rather than from an individual parton". In other 
words the distributions of the produced hadrons "retains the memory" not only of the 
struck quark type but also of the target remnant and, hence, the entire string configura- 
tion. 

2 The option LST(8)=1 used in reference [64] of ^ corresponds to inclusion in the generation of the 
first order QCD matrix elements for gluon radiation and photon-gluon fusion. With this option Eq. Q 
has to be replaced with the NLO expression to include the gluon distribution and fragmentation functions 
which arc not involved in the LO purity analysis. However, it was checked that the conclusions of the 
present paper remain valid independently of the choices of the LEPTO settings at HERMES energy. 
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Figure 1: "Quark fragmentation functions" to tt + calculated for HERMES conditions. 
Solid line - with cut x < 0.094, dashed line - with cut x > 0.094. Fragmentation functions 
from Kretzer: dotted line -forQ 2 = 1.5(GeV/c) 2 , dot- dashed line -forQ 2 = 3A(GeV/c) 2 . 



Event generation in LEPTO proceeds in three steps: first the hard scattering kinematics 
(x, Q 2 ) is chosen from the differential DIS cross-section. Second, the struck quark flavor is 
chosen. Third, the string is set up and hadronized according to LUND model implemented 
in JETSET program [3] . Within this approach the SIDIS cross section can be expressed as 

a h N {x, z, Q 2 ) oc (1 + (1 - y) 2 ) £ e\ q(x, Q 2 ) H% /N {x, z, Q 2 ) , (6) 

where the functions H^ N (x, xf, Q 2 ) are describing the conditional probability of hadron h 
production in the hadronization of the system formed by struck quark q and corresponding 
target remnant. Let us call them hadronization functions, HFs. The target remnant type 
and, hence, the whole fragmenting system configuration depends both on the nucleon 
type and on the struck quark type 3 . This means that in contrast with independent 

3 The fragmenting system in LEPTO/ JETSET is in general more complicated than a quark-diquark string. 
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HERMES 




fragmentation functions the HFs are non universal - they depend on the process type and 
energy. 

Note, that Eq. © is valid not only in the current fragmentation region but in the 
whole xf interval. 

The product q(x, Q 2 ) H\ N (x, z, Q 2 ) is the probability to find the parton q in the nu- 
cleon, N, and, after hard interaction, to create a hadron h in the string hadronization. By 
its physical meaning (probabilistic interpretation) this object represents nothing else but 
the fracture functions discussed in There exist certain arguments based on handbag 
diagram dominance that this concept may be applied even in the current fragmentation 
region of SIDIS 6 . The LEPTO/JETSET Monte Carlo program can be considered as a 
model for these functions. It is clear that from the generated samples one can actually 
extract only HFs and as one can see from Fig.^and Fig. Eleven in the current fragmenta- 
tion region one cannot neglect the dependence of this functions on the Bjorken x variable 

The target remnant state depends on the removed active parton type and the whole fragmenting system 
may contain multi-string configurations |2J [3] . 
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and on the target type. 



HERMES 
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Figure 3: "Quark fragmentation functions" dependence on diquark type calculated for 
HERMES conditions with cut x > 0.1. Solid line - scalar diquark, dashed line - vector 
diquark. 

As a further consideration of the target remnant state influence on the "quark frag- 
mentation functions" let us now consider the spin inside LEPTO. In the simplest case 
when the valence w-quark is removed by hard scattering from a proton, the target remnant 
is a scalar, (ud)o, diquark with probability wq = 0.75 or a vector, (ud)i, diquark with 
relative probability W\ = 1 — wq = 0.25 [21 Ej- In Fig. El the "fragmentation functions" for 
h = 7r + , 7r~, K + and K~ extracted with x > 0.1 cut are presented for the cases when the 
target remnant diquark is chosen to be a 100% scalar (w = 1) or a vector one (wo = 0). 
We see that HFs calculated with our generated samples exhibit dependence (at 5-10% 
level) on the target remnant spin state already in unpolarized SIDIS. 

This dependence has a quite general origin. To understand the underlying physics let 
us consider, as an example, a K + meson production. At some stage of hadronization the 
ss-pair is created and s-quark can combine with the remnant diquark to form a strange 
baryon. The formation of first rank A hyperon is possible for a scalar diquark remnant 
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and forbidden for the case of vector diquark. In the second case only a heavier strange 
hyperon can be formed at first rank. Then, due to the energy-momentum conservation 
the available energy for a K + meson production will be less in the second case compared 
with the first case. This is the reason why in Fig. |3]the "w-quark fragmentation functions 
into K +n is higher for scalar diquark sample. 

3 Polarized SIDIS and String Fagmentation 

In the ordinary factorized picture the same unpolarized fragmentation function is entering 
into expression for the cross section of unpolarized and polarized SIDIS. Let us now 
compare that what will happen if we include spin and generalize the picture with the HFs 
for SIDIS. 

At present there is no the polarized version of the string fragmentation Monte-Carlo 
program for event generation. It is evident that the description of polarized SIDIS is 
more complicated than in the unpolarized case. Let us, as an example, again consider 
the simplest case when the valence w-quark with positive or negative (u~) helicity 
is removed from nucleon with positive or negative helicity, N + or N~ . Within the SU(6) 
quark-diquark model used in Ej the polarized nucleon wave functions are given as 
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n 



-.{u + [3(ud) Q ,o + (W)i,o] - y/2u~(vd)i tl - V2d + (uu) lfi + 2d~(uu) 1)1 }, (7) 



-={d + [3{ud) 0fi + [ud) lfi ] - V2d-(ud) 1}1 - V2u + (dd) 1>0 + 2u~(dd) ltl }, (8) 
/18 



where (qiq2){i,k) stands for the diquark formed by the q±- and g 2 -quarks with spin i and 
spin projection k. 

Using the explicit form of the polarized nucleon wave functions one can calculate the 
relative probabilities, w, of different target remnant states (and hence the states of entire 
string) depending on the struck quark and the nucleon polarizations. For example, when 
the M + -quark is removed from the p + we get the following string configurations with 
corresponding probabilities w 

/e^f {HVo u+h w = - 9 ' (9) 

[{(ud) 1>0 u + }, w = 0.1, 

where {{qiq2)i,k <? + } denotes the string formed by the struck quark q + and the diquark 
(qiq2)i,j- Similarly, when the u + -quark is removed from p~ we get 

p- Qu + ^> {(ud)^ u + }, w = l. (10) 

For the neutron target we have 

n + u + {{dd) lfi u + }, w = l, (11) 

and 

rr eu + {(dd)i,-i u + }, w = l. (12) 
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The relations (Pll2j) demonstrate that the string configuration indeed depends not only 
on the struck quark type and its polarization but also on the target type and polarization. 
As we have seen in Sec. |2]the description of hadron production in the current fragmentation 
region of SIDIS within the LUND fragmentation model does not correspond exactly to the 
commonly adopted simple picture of independent quark fragmentation but rather to the 
more complicated approach based on fracture functions. Even in unpolarized SIDIS HFs 
depend on target fragment spin states as it is demonstrated in Fig. Then in polarized 
SIDIS the dependence on the target and struck quark polarizations appears. So, one has 
to generalize Eq. (JHJ) for the polarized SIDIS case. 

Let us start with the SIDIS cross section ct NXiXn for the positive helicity lepton, A; = +1 
and hadron, \ N = +1: 

a h N++ cx H h q/N++ + (1 - yfq- H h q/N _ + }, (13) 

i 

where H\ NX Xn describes the production probability of the hadron h in the quark-target 
remnant system fragmentation and depends not only on x and z (or x and xp) but also 
on the struck quark and nucleon helicities, X q = ±1 and = ±1. Similarly 

a h N+ _ oc J24{ ( 1~Kn + - + (1 " y?q + H h q/N _„}. (14) 
g 

The partially polarized beam state, l Xl , (with helicity A;) can be described as l Xl = 
1/2(1 + Xi)l + + 1/2(1 — A/)Z~ and similarly for the nucleon. Then in the general case of 
arbitrary polarized beam and target we have 

<x t x N oc (1 + A,)(l + X N )a h N++ + (1 + A,)(l - X N )a h N+ _+ 

(1 - A,)(l + Ajv)<4_ + + (1 - A,)(l - X N )a h N __. (15) 

Now, using the relations H^ N++ = H^ N and H^ N+ _ = H^ N + which follow from 

parity invariance and introducing 

zrh Tjh , zrh 

n q/N — n q/N++ n q/N+- 

A-ffg/Ar = H q / N++ — H q/N+ _ (16) 

after simple algebra one get: 

[1 + (1 - yf] e 2 q {qH h q/N + ^H h q/N }+ 
i 

X { X N \\ - (1 - yf] J2e 2 q {AqH^ /N + qAH^}, (17) 
q 

where now A; and Xn are the (arbitrary) beam and target helicities. 

Eq. (JT7j) is very similar to the equation proposed in [7] . The difference is that functions 
Hffi N and AH^ N are not independent quark fragmentation functions like in [7j. It is well 
known that the single spin dependence is forbidden by the parity invariance for integrated 
over transverse momentum independent fragmentation. By this reason the same quark 
fragmentation functions enter in Eqs. (^EJ) .In contrast, HFs describe the probability of 
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hadron production in the hadronization of the whole struck quark - target remnant system 
and, hence, one deal with the double spin effects. Even for hadrons produced in the current 
fragmentation region these functions can depend not only on the fraction of the quark 
energy carried by produced hadron but also on the whole hadronic CMS energy and the 
target and the struck quark polarizations. 

When integrated over the whole available phase space of the selected hadron, Eq. (JT7jl 
transforms to the standard parton model expression for the polarized DIS, provided that 
the following sum rules hold valid: 



J2 h JdzzH* N (x,z,Q 2 ) = l, 

J2jdzzAH^ N (x,z,Q 2 )=0. (18) 

4 Discussion and Conclusions 

The standard expression for the SIDIS description in the current fragmentation region 4 
is obtained if one assume that 

H; /N (x, z, Q 2 ) - D h q (z, Q 2 ) (19) 

and 

AH^ /N (x,z,Q 2 )^0. (20) 

As we have demonstrated in Sec. El relation (|19|) is not correct for the HERMES 
experimental conditions in the LUND fragmentation model. On the other hand we have 
seen in Fig. |3] that hadronization in this model depends on the target remnant spin 
quantum numbers. In the case of polarized SIDIS the relative probabilities of different 
target remnant states depend on the target and struck quark polarizations, see Eqs. (0- 
ITU]) . So, there is no any reason to believe that relation ()20|) will hold for the polarized 
SIDIS at moderate energies. Thus, a new nonperturbative inputs - the polarized HFs, 
AHy, N (x, z, Q 2 ), are needed. In this case Eq. (jSJ) underlying the purity method is not 
exact and must to be replaced by 

V e 2 a(r 2 )H h (r z O 2 ) \ Aq(x ' Q2) I AH *M X > Z > Q2) ) 

Ai{X,Z,Q)= Q 2 )AH h (x z O 2 ) ' 

E a 4 „(», q>) flj w (,, Q*)[i + "^^Al ] 

with extra contributions in the numerator and denominator as compared to Eq. (j3J). As 
one can see from Fig. |3] the second term in the square brackets in numerator can reach 
5-10%. Its influence can still be negligible in the unpolarized cross section (second term 
in square brackets of denominator), since it entering multiplied to quark polarization. 
Whereas in numerator it can be comparable with (or even greater than) the first term. 

Though we are not able to calculate from the first principles the HFs, one can try 
to estimate the possible effects of polarized HFs on the extraction of polarized quark 

4 Note that these expressions are based on QCD factorization theorems, which represent asymptotic 
statements valid for very high lepton energies, Q 2 and W. Only under these conditions can one neglect 
the influence of the target remnant on hadron production in the current fragmentation region. 
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Figure 4: "Quark fragmentation functions" to 7r + calculated for COMPASS conditions. 
Solid line - with the x < 0.028, dashed line - with cut x > 0.028. Fragmentation functions 
from Kretzer: dotted line -forQ 2 = 1.7(GeV/c) 2 , dot- dashed line -forQ 2 = 5.0(GeV/c) 2 . 

distribution. For example, one could estimate the effects of scalar and vector diquarks 
using formalism of Sec. 121 and calculate TT n /N } ' fi? for each x-bin. 

In a model for the extra contribution in the numerator of Eq. (}2~T]) has been consid- 
ered and it was demonstrated that our ignorance of polarized HFs may lead to incorrect 
results for polarized sea quark distributions. Here it is demonstrated that using LEPTO in 
analysis based on the factorized approach is inconsistent. As a consequence it is essential 
in using factorization approach to polarized quark distribution extraction first to check it, 
since the polarization dependence is more sensitive to factorization that the unpolarized 
multiplicity distributions. 

It is interesting to study how these effects depend on the energy and in particular how 
they will affect the COMPASS [8 a analysis. In Fig. 0] the same distributions as in Fig. Q 
are presented for COMPASS kinematics. One can see that the dependence on Bjorken 
variable of the "quark fragmentation functions" extracted from generated samples is less 
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pronounced that in case of HERMES experimental conditions. The same observation is 
also valid for the dependence upon the target type and the target remnant spin state 
dependencies. This means that that polarized HFs may be negligible in the current 
fragmentation region at high energies. 

The string configurations considered in ([DlfP2|) correspond to the simplest case of re- 
moving the valence quark from nucleon. In the case, when the virtual photon interacts 
with the sea quark or higher order hard scattering processes are considered, the target 
remnant and the final parton configuration are more complicated [2|. For example, in 
photon-gluon fusion the target remnant is split into a quark and a diquark that form 
two respective separate strings with the antiquark and quark produced in the fusion pro- 
cess. One can generalize Eq. (|17|) to include the corrections from higher order QCD hard 
processes. The generalized expression for the polarized cross section of single and two 
hadron productions will, in addition, contain new unknown HFs, AH^ N (x, z,p!^, Q 2 ) and 

AHgj£ 2 (x, Zi, Z2,Py,Py 2 , Q 2 ), with the corresponding distribution functions, (A)g(x, Q 2 ). 
This means that the validity of the Monte Carlo based approach JUj to extract the po- 
larized gluon distribution may also be questionable at moderate energies. 

It has been recently noted that the appearance of separate distribution and frag- 
mentation functions cannot be proven in general, but is rather assumed and justified 
a posteriori, while the natural framework to describe SIDIS involves fracture functions. 
These functions can be also generalized to describe the T-odd single spin asymmetries jH] . 
As it was mentioned in Sec. 121 within the LUND fragmentation framework the fracture 
functions can be represented by products of distribution functions and HFs. Recent de- 
velopments in the theory of SIDIS for single spin asymmetries and diffractive phenomena 
also show that one cannot neglect the interaction of (colored) removed partons and target 
remnants (see, for example, JT] and references therein, and [S])- This again indicates 
that the description of SIDIS based on the naive parton model with the independent 
fragmentation is only an approximation, to be justified at moderate energies. 

Let us stress that the approach developed in Sec. El is not directly derived from QCD 
but is similar to one used in LEPTO event generator. It assumes that in the DIS regime 
the hard scattering and the hadronization are factorized. As it was mentioned long time 
ago [12 the concept of independent fragmentation can be justified only when there is 
enough phase space for the final hadronic system. The important issue is not only high 
Q 2 but also large rapidity interval available for a given hadron production. In contrast 
to the independent fragmentation model the LUND model deals with the whole final 
quark - target remnant hadronization and takes into account the energy-momentum 
conservation, the color flow and the quantum number correlations. As a consequence, 
at moderate W there is a non-negligible influence of the target remnant state on the 
distributions of hadrons produced in the current fragmentation region, as it explained in 
the end of Sec. 121 
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